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ABSTRACT 

Aims. Our goal is to understand the nature of blazars and the mechanisms for the generation of high-energy y-rays, through the 
investigation of the blazar 3C 66A. 

Methods. We model the high energy spectrum of 3C 66A, which has been observed recently with the Fermi-LAT and VERITAS 
telescope. The spectrum has a hard change from the energy range of 0.2-100 GeV to 200-500 GeV in recent almost contemporaneous 
observations of two telescopes. 

Results. The de-absorbed VERITAS spectrum greatly depends on the redshift, which is highly uncertain. If z=0.444 is adopted, we 
are able to use the SSC model to produce the Fermi-LAT component and the EC model to the VERITAS component. However, if 
z=0.1, the intrinsic VERITAS spectrum will be softer, there will be a smooth link between the Fermi-LAT and VERITAS spectra 
which can be explained using a SSC model. 
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1. INTRODUCTION 

Blazars are a peculiar class of active galactic nuclei (AGN) and 
their jets point at small angles with respect to our line of sight. 
Many of them have been observed at all wavelengths, from ra- 
dio to very high energy (VHE) y-rays. Their spectral energy dis- 
tribution (SED) consists of two bumps which are attributed to 
the synchrotron and the inverse Compton (IC) emission of ultra- 
relativistic particles. The different soft photon sources deduce 
synchrotron self-Compton (SSC) and external Compton (EC) 
models to produce high energy emission. In the SSC model, the 
soft photons are provided by the synchrotron emission of the 
■ same electrons dMarscher 1980I IGhisellini & Maraschi 1989I 
Mars cher & Travis 1996l l. However in the EC model, the soft 
photons mostly come from the outside of the jet, such as outer 
disk, broad-line region (BLR) clouds (etc. IDermeret al. 19931 

ISikora et al. 1994b . 

In 3C 66A, |Miller et al. (1978)| have given the redshift 
z=0.444 by a weak Mg II emission line detection, but it is 
very uncertain (Br amel et al. 2005T l. When 3C 66A locates at 
z=0.444, its TeV photons will suffer the strong pair produc- 
tion absorption of the Extragalactic Background light (EBL). 
After corrected by the EBL absorption, TeV emission presents 
an inverted intrinsic spectrum (see the Acciari et al. (2009) Fig. 

2, the de-absorbed photon spectral index is calculated to be 
1.1 ± 0.4). In this paper, we take z = 0.444 to analyze TeV 
emissive mechanism, and will discuss the behavior of the de- 
absorbed VERITAS spectrum in different redshifts. Generally, 
3C 66A is classified as a low-frequency peaked BL Lac ob- 
ject (LBL). The peak of the low-frequency component of LBLs 
usually lies in the IR or optical regime, whereas the peak 
of high-energy component locates at several GeV. The lumi- 
nosity of y-ray is typically comparable to or slightly higher 
than the luminosity of the synchrotron component. Such as, 



|Joshi & Bottcher (2007)| have argued that for 3C 66A the peak of 
low-frequency component locates at the optical regime, the peak 
of high-frequency component reaches multi MeV-GeV range. 
However, Petri et al. (2003) have revealed that the synchrotron 



peak locates in between 10 13 and 10 16 Hz, then 3C 66A is clas- 
sified as an intermediate-frequency peaked BL Lac (IBL). From 
the X-ray spectrum with the photon spectral index F ~ 2.5 
(Bottcher et al. 2005; Donato et al. 2005; Foschini et al. 2006), 
which might be the tail of the synchrotron emission, 3C 66A 
is considered as an IBL in this paper. 3C 66A is observed in 
radio, IR, optical, X-rays and y-rays, and shows large luminos- 
ity variations. As described in Bottch eFet al. (2005)| the object 
exhibits several outbursts in the optical band and the variations 
of Am ~ 0.3-0.5 over a timescale of several days. Until now, 
the majority of BL Lacs detected at VHE (very high energy: 
E > 100 GeV) are HBLs (high-frequency peaked BL Lacs). 
Only IBL W Com ae dAcciari et al. 20081), LBL B L Lacertae 
dAlbert et al. 20071 > and 3C 279 dAlbert et al. 20081) display the 
potential to enlarge the extragalactic TeV source. For 3C 66A, 
the Crimean Astrophysical Observatory have reported a 5. 1 cr de- 
tection above 900 GeV(Stepanyan et al. 2002). 

Recently, VERITAS have carried out 14 hours' observa- 
tions for 3C 66A from September 2007 through January 2008 
(hereafter, the 2007-2008 season), and from September through 
November 2008 (hereafter, the 2008-2009 season) a further 46 
hours' data have been taken dAcciari et al. 20091 . Due to the 
limited spatial resolution of Cherenkov telescopes it is diffi- 
cult to accurately identify the emission region. The radio galaxy 
3C 66B lies in the same view field of 3C 66A at a separa- 
tion of 0.12° and is also a plausible source of VHE radiation 
dTavecchio & Ghisellini 20081 ). The recent detection by MAGIC 
favored 3C 66B as VHE source and excluded 3C 66 A at an 85% 
confidence level dAliu et al. 20 09). However , VERITAS have 
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found that 3C 66A lies 0.01° from the fit position while 3C 66B 
lies 0.13° away, and 3C 66A is VHE source. If 3C 66A has a 
redshift of z = 0.444, its de-absorbed spectral index is 1.1 ±0.4 
showing very hard intrinsic spectrum dAcciari et al. 20091 . In 
first three months, the Fermi-LAT Gamma-ray Space Telescope 
have observed 3C 66 A (Abdoetal. 2009), almost at contem- 
poraneous observation with VERITAS in the 2008-2009 sea- 
son. However, very soft spectrum with the spectral index of 
1.97+0.04 appears in the Fermi-LAT observing energy range. 
The y-ray spectrum suddenly hardens from 0.2-100 GeV to 200- 
500 GeV and challenges the one-zone homogeneous SSC model. 

In Section 2 we present the jet models for application to 
3C 66A. We use the observed data to constrain the model pa- 
rameters in Section 3. We finish with discussions and conclu- 
sions in Section 4. Throughout this paper, we use a soft cosmol- 
ogy with a deceleration factor qo = 0.5 and a Hubble constant 
H = 15kms- y Mpc- 1 . 



2. THE MODELS 

We use a public model of Georganopoulos ~et al. (2007)| to de- 
scribe the observed spectrum of 3C 66A. The emission region is 
assumed to be a sphere (blob) with radius R, permeated by a ho- 
mogeneous magnetic field B. The blob moves with bulk Lorentz 
factor F through an external photon field with a black body spec- 
trum, at an angle 9 (in this work 8 ~ 1 /T is assumed) with re- 
spect to the line of sight, and has a Doppler factor 5 ~ T. The 
relativistic electrons are continuously injected into the blob at 
a rate Q e (y) = Qoy~ s [ernes' 1 ] between y min and y max , where 
y m i n is the minimum Lorentz factor of the injected electrons and 
should not be confused with the minimum Lorentz factor of the 
emitting particles. The injection correspond to a luminosity L; n j, 



|Geor ganopoulos et al. (2007) modified an analytical appro xima 
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Where V = 4nR 3 /3 is the volume of the blob. Injected elec- 
trons might obtain energy from acceleration processes before 
into the emission region, such as diffusive shock acceleration, 
second order Fermi acceleration and gradual shear acceleration. 
Detailed discussions about the acceleration processes are beyond 
the scope of this paper. The electrons suffer synchrotron and in- 
verse Compton losses and eventually escape from the emission 
region. 

The time-dependent evolution of the electron population 
N(y, t) [cm -3 ] inside the emission region is governed by, 



tion used by |Moderski et al. (2005)| y ic = ^y 2 U I F KN , where 



{/, is the energy density of seed photons including synchrotron 
photons and external photons (such as the reprocessed pho- 
tons by B road Line Region (BLR ) dSikora et al. 19941 l). F KN is 
given by Moderski et al. (2005)[ In this paper, the latter pho- 
tons are assumed to be a blackbody radiation with peak fre- 
quency v ext and energy density U ext (all seen in the observer 
frame). For the beaming of the EC emission, we use the recipe 



of |Georganopoulos et al. (2001) 



3. MODELING PARAMETERS 

First of all, we use the observed quantities to estimate the 
physical parameters in the blob, and then use these values to 
reproduce the observed SED. [Mastichiadi s & Kirk (1997)1 an d 
|Konopelko et a l. (2003) have estimated the parameter relations 
of the inverse Compton scattering in Klein-Nishina regime in a 
homogeneous SSC scenario. However, Paggi et al. (2009) have 
found that the simple relations of parameters in Thomson regime 
satisfy the observations of LBL or IBL sources, and that the re- 
lations in the extreme KN limit is not suited to HBLs, indicating 
that the inverse Compton scattering for HBLs just borders the 
KN regime. For 3C 66A, the Fermi-LAT spectra showing flat and 
rising shapes and higher energy observations of VERITAS both 
indicate that the IC scattering of Fermi-LAT spectra do unlikely 
enter into the Klein-Nishina regime. Therefore, we assume that 
the Fermi-LAT spectrum just cover the peak of SSC emission, 
and use the parameter relations in Thomson regime to rudely 
estimate the parameters. In fact, we consider the Klein-Nishina 
effect in producing the high energy spectra. 

The size of the emitting region, R, can be constrained by 
measurements of variability timesscales, i.e. R< 6ct var /(l + z). 
The observations of VERITAS have shown the variability to 
be the time-scale of days ( Accia net al. 20091 . The multiwave 
campaign of 3C 66A by the Whole Earth Blazar Telescope 
(WEBT) in 2007-2008 have observed several bright flares on 
time scales of ~ 10 days (Bottcher et al. 2005). These observa- 
tions suggest the size of the emitting region to be about 10 16 cm 
if we assume 6 is order of 10. From the peak frequencies of 
synchrotron and inverse Compton radiation, we can estimate 
the y pea k (i.e. the Lorentz factor of the electrons emitting at 
the peaks of the synchrotron and SSC components) follow as 
(ITavecchio et al. 1998k 
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Here, (dy/dt)^ is the radiative energy loss rate, due to syn- 
chrotron, SSC and/or EC emission. f esc is the electron escape 
timescale which is several light crossing times, in this work 5 
R/c is adopted. 

The energy loss rates of electrons caused by the syn- 
chrotron are given by: y s = j^.7 2 Ub, where Ue - B 2 /(8n) 
is the magnetic energy density. The energy loss rates of in- 
verse Compton emission, fie, have excellent analytical ex- 
pressions for the Thomson regime x <k 1 and for the deep 
Klein-Nishina (KN) regime x » 1, but have not ones for the 
middle regime of i = ey (e is the energy of the incoming 
photon in units of the electron rest mass). To overcome this, 



and we obtain 



and 
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where, vi = eBI2mn e c =* 2.8 x 10 6 B Hz. From the equation 

L SSC L SSC tpeak,SSC U' s L s 8tt 2L s 
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we have 
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where L s and L$sc are me observed total luminosity of the 
synchrotron peak and SSC peak, L ssc and L s are the lumi- 
nosity in the comoving frame, and y pea k,s and y pea k,ssc denote 
the synchrotron and SSC cooling rates of the electrons, U s and 
U B are the comoving energy densities of synchrotron photons 
and magnetic fields, the di is the luminosity distance, F vs and 
F v ,ssc are the energy flux at peaks. Using the equation (O and 
(O and R « dct var /(l + z) (here, t var ~ 2 days is assumed), 
we are able to derive the 6 and B from the v s , vssc* v s F VtS , 
and v S scF v ,ssc- We take v^4x 10 15 Hz (Perri et al. 2003), 
VsFv,s ~ 7 x 10~ n erg ■ cmT 2 s~ l estimated from historical data, 
v ssc ~ 1-21 x 10 25 Hz which is the middle band of Fermi-LAT, 
vsscF v ,ssc ~ 4.9 x 10~ n erg ■ cmr 2 s~ l which is an approxima- 
tive value observed by the Fermi-LAT, and z=0.444. We then 
get 5 ~ 29 and B~ 0.025G. 6 can also be estimated by other 
methods. The apparent speed B app -\2.\c has been presented by 
Bottc her et al. (2005)| Modeling the non-simultaneous SED of 
3C 66A, Ghisellini et a l. (1998) suggested T -14, which is a 
typical value of blazars. Jor stad et al. (200 1)| favored a high su- 
perluminal speed up to B app » 27c and got F > 27 or 6 ~ 30. 
However, it is noted that we use the flux of the Fermi-LAT 
band to estimate the parameters, if the modeling spectrum in- 
cludes the VERITAS spectrum, the estimated parameters will be 
modified. In fact, modeling quasi-simultaneous observations of 
Fermi-LAT and VERITAS obtains 5 = 35 (see the Table 1.). 
The steady state distribution of emitting electrons is given 

by: 

N(y) = * : . (9) 

r 

What is the criterion for steady state? The code of 



Georganopoulos et al. (2007) firstly takes a time step equal to 
the cooling time of the lowest energy electrons under the syn- 
chrotron and external Compton losses, and then calculates N(y) 
to a steady state where the injected electrons have fully been 
cooled, through an adaptive time step. 

We emphasize the selection of the spectral index s of injected 
electrons. Based on the X-ray photon spectral index of T * 2.5 
(Bottcher et al. 2005; IDonato et al. 20051 IFoschini et al. 2006D , 
we deduce the spectral index of emitting electrons to be 4. 
Assuming the observed X-rays to be from synchrotron emis- 
sion of the cooled electrons, the spectral index of the injected 
electrons is given by s = 3. |Celotti & Ghisellini (2008)1 adopted 
s — 3.6 to model the SED of 3C 66A. Therefore, we use s ~ 3 to 
model the observed data. The particle injection, radiative cool- 
ing, and escape from the emission region might yields a tempo- 
rary quasi-equilibrium state described by a broken power-law. 
The balance between escape and radiative cooling will lead to a 
break in the equilibrium particle distribution at a break Lorentz 
factor jb, where t esc = t coo i(yi,). The cooling time scale is eval- 
uated taking into account synchrotron, SSC and EC cooling. 
Depending on whether y/, is greater than or less than y, m „, the 
system will be in the slow cooling or fast cooling regime. In the 
fast cooling regime (y^, < y m ,„), the equilibrium distribution will 
be a broken power-law with N(y) oc y~ 2 for yt, < y < y„,i„ 
and N(y) oc for y mi „ < y < y mnx . In the slow cool- 

ing regime (yi, > y mlM ), the equilibrium distribution will be 
N(y) oc y- 2 for y < y min , N(y) oc y- 5 for y min < y < y b 
and N(y) oc for yi, < y < y max . Since for thin syn- 

chrotron emission the energy spectral index is related to that of 



the emitting electrons as a = (q - l)/2 , where q is the spectral 
index of the emitting electrons distribution, and it is immedi- 
ate to see that the peak in the vF Y spectrum occurs where q=3. 
When s > 3 (in this work s=3.4 and 3.3 are adopted), y,„,„ (the 
minimum Lorentz factor of the injected electrons) almost cor- 
responds to y pea k dGhisellini et al. 1998 ; Ghisell ini et al. 2 002). 
Jmax presents the balance between the acceleration and cooling 
and has small impact upon the SED, it is usually taken to be 
10 5 - 10 7 (Inoue & Takahara 1996). 

Considering the EC emission of electrons, we need to esti- 
mate U ext and v ext by mainly considering the soft photons re- 
processed by the BLR. v ext is usually considered to be around 
optical-UV wavebands, we let v ex , = 2.5 x 10 15 Hz. The 
L ext of FSRQ is easily estimated from emission line spec- 
tra or UV-excesses, while it is difficult to estimate the L ext 
of BL Lacs. Assuming the luminosity of accretion disk to be 
Lq = 1 x 10 45 erg ■ s , which is larger than the jet lumi- 
nosity and does not produce a blue bump in the simulated 
SED (Joshi & Bottcher 200 7), and taking R ext ~ 10 17 cm 
(ITavecchio & Ghisellini 20081 . we get the upper-limit U exttU = 
2.65 x 10~ 2 erg • cm~ 3 assuming the reprocessing efficiency of 
the BLR to be 0.1. In this work, we estimate U ex t through mod- 
eling the spectrum under the condition of U ex , < U exlM . Taking 
U ext = 3.5 x 10~ 6 erg ■ cm~ 3 , we can already reproduce the ob- 
served data of VERITAS. Therefore the soft photons reprocessed 
by the BLR can provide the EC emission of electrons to produce 
high energy radiation. 

In the Fig. [1] we present the modeling results for 
the observed data of 3C 66A. The solid triangles are 
the data of |Perri et al. (2003)] the open circles come 

blue bow-tie show 
squares denote the 



from [Bottcher et al. (2009) 
the Fermi-LAT data, and 



the 

the green 

VERITAS results corrected by EBL absorption according 
to [Franceschi ni et al. (2008)| In the figure, we include a 
quasi-simultaneous data including the near-infrared, optical, 
UV-Optical, and X-ray observations (Reyes et al. 2009). We use 
the synchrotron emission model to model the lower energy part 
(black solid line), and use the SSC model to reproduce the spec- 
trum observed by the Fermi-LAT (red dash line). Particularly 
we use the EC model to model the VERITAS spectrum (green 
dot line). It is indicated that the harden spectrum from the 
Fermi-LAT to VERITAS energy ranges could exhibit an EC 
spectrum. 

4. DISCUSSIONS AND CONCLUSIONS 

The redshift of 3C 66A has an uncertain value 
(Bra mel et al. 2 005). and is usually adopted as z = 0.444. 
However the redshift is crucial in constructing intrin- 
sic high energy spectrum because of the EBL absorption 
(IHauser & Dwek 200 H . This absorption decreases the observed 
flux and softens the observed spectrum. If the redshift is less than 
0.444, such as just > 0.096 suggested by |Finke etal. (2008)| the 
intrinsic spectrum in the VERITAS energy range will be softer. 
There will be a smooth link spectrum between the Fermi-LAT 
and the VERITAS energy ranges. In the Fig. [2] we generate 
the intrinsic spectra of VERITAS observations corrected by 
EBL absorption according to Franceschini et al. (2008) model, 
assuming the source to be at the different redshifts z=0.03, 0.1, 
0.3, and 0.5. It is shown that the de-absorbed spectrum strongly 
depends on the redshift. When z=0.3, the de-absorbed spectrum 
has a little inverted, but it becomes an inverted spectrum in 
z=0.5. If the redshift is less than 0.1, the de-absorbed spectrum 
will present the usual SED of an LBL such as W Comae 
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Table 1. Parameters for the SSC+EC or SSC model used to re- 
produce the SED (Fig.Q]and Fig. [3). 



Parameters 


If z=0.444 


If z=0.1 


L mj (10 41 ergs-') 


4.0 


6.3 


Tmin 


2.5 x 10 4 


1.5 x 10 4 




8.0 x 10 5 


8.0 x 10 5 


S 


3.4 


3.3 


R (10 16 cm) 


7.0 


7.0 


6 


35 


30 


B (Gauss) 


0.034 


0.05 


Vext (HZ) 


2.5 x 10 15 




Uext (erg/cm 3 ) 


3.5 x 10- 6 





dAcciari et al. 20081 . In the Fig. [3] we show the de-absorbed 
SEDs under z=0.1 and the modeling. A smooth spectrum 
can link the Fermi-LAT and VERITAS data and be explained 
with a SSC model, in which U B = 9.95 x \Q~ s erg ■ cm' 3 and 
U e =2.31 x IQ- 3 erg -cm.- 3 . 

In fact the bulk motion of the emitting blob affects the ob- 
served SED (e.g., Derm er 19951 Georganopoulos et al. (2001)). 
The peak frequencies are given by v, ~ 3.7 x lQ 6 y 2 peak 6B/(l +z) , 



vssc ~ \y 2 peak v s ,wdv EC 



\y 2 peak V ^ T6 l^ We Can See that 



Vec would be larger than v ssc provided the blob has larger 6 or 
T (Usually the viewing angle 6 ~ 1/r is assumed, 6 ~ F.) From 

the ratio of peak luminosity, 7^- * ^F- = i^S. « ^ « 

L J LsSC L ssc jpeak.SSC U^„ 

^4p2f^ff^ (U ext is amplified by a factor of F 2 , see the Sikoraet 
al. (1994) and Dermer (1995)), where £ is the reproduce fraction 
of the L ext , we show that j^- is strongly affected by the bulk mo- 
tion of the blob. In the Fig. 1, using U B = 4.6 x \0~ 5 erg ■ cm~ 3 , 
U' e = 1.67 x \Q~ 3 erg ■ cm' 3 , and U' ext = 3.5 x 10~ 6 erg ■ cmT 3 we 
can model the SED. U ext is obviously lower than U B , however, 
the EC luminosity is comparable with the synchrotron ones (see 
the Fig. 1). In the observer frame, the beaming factor is differ- 
ent for EC (<5 4+2a dDermer 19951 Georganopoulos et al. 2001)), 
synchrotron and SSC emission (6 3+a ). The difference of EC and 
synchrotron luminosity is reasonable by considering their beam- 
ing factor. 

The EC emission is less clear for the BL Lac objects. The 
lack of strong emission lines and UV excesses suggests that the 
external photon density %L ext is very low, while the Lorentz fac- 
tor of BL Lac objects is typically smaller than that of quasars 
dPiner et al. 20081 1 . Their high energy emission strongly favors 
the SSC mechanism over the EC mechanism. But, 3C 66A might 
be an exception and reveal an existence of larger bulk velocity 
in the high energy emissive region. Therefore, the high energy 
emission caused by EC mechanism is likely observed in the IBL. 
This possibility needs future Fermi-LAT and VERTAS observa- 
tions for 3C 66A and a precise redshift determination. 
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